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QQ ■ We study the evolution of a group of TRACE 195 A coronal loops overlying a re- 

verse S-shaped filament on 2001 June 15. These loops were initially pushed upward with 

. the filament ascending and kinking slowly, but as soon as the filament rose explosively, 

^/^ . they began to contract at a speed of ~100 km s~^, and sustained for at least 12 min, 

^ i presumably due to the reduced magnetic pressure underneath with the filament escap- 
ing. Despite the contraction following the expansion, the loops of interest remained 

O I largely intact during the filament eruption, rather than formed via reconnection. These 

' contracting loops naturally formed a shrinking trap, in which hot electrons of several 

, keV, in an order of magnitude estimation, can be accelerated to nonthermal energies. 

A single hard X-ray burst, with no corresponding rise in GOES soft X-ray flux, was 

^ , recorded by the Hard X-ray Telescope (HXT) on board Yohkoh, when the contracting 

! loops expectedly approached the post-flare arcade originating from the filament erup- 

\ tion. HXT images reveal a coronal source distinctly above the top of the soft X-ray 

^—^ I arcade by '^15". The injecting electron population for the coronal source (thin target) 

OO ' is hardening by --^1.5 powers relative to the footpoint emission (thick target), which is 

' consistent with electron trapping in the weak diffusion limit. Although we can not rule 

■ out additional reconnection, observational evidences suggest that the shrinking coronal 

^ . trap may play a significant role in the observed nonthermal hard X-ray emission during 

^ , the flare decay phase. 

Subject headings: Sun: flares — Sun: X-rays, gamma rays — Sun: filaments 



1. Introduction 



In the standard flar e model (CSHKP model: ICarmichaellll964l : ISturrocklll966l : lHiravamalll974l : 
Kopp &: Pneumanlll976l ). the initial driver of the flare process is a rising filament above the polarity 
inversion line. The rising filament stretches the overlying field lines, resulting in the formation of 
a current sheet underneath and the consequent cancelation of magnetic fluxes of opposite polarity. 
It is energetically impossible, however, to open up the closed magnetic fiel d lines via a purely 
ideal-MHD process in the force-free configuration (the Aly-Sturrock paradox: lAlylll984l : ISturrock 
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199lh . One way to avoid th is constraint is to open only a portion of the closed field lines (see 
the review by Lin et al.1 120031 ) . Considering a flux rope confined by an overlying magnetic arcade, 
Sturrock et al] ~( 2001 ) demonstrated that it is energetically favorable for part of the rope to erupt 
into interplanetary space. 



This partial eruption concept has been further inv estigated both observationall y (e.g., 



2000 . 12001 
2004 



Fanlboosi : Icibson k Fan! 



2007l:lLiu et al.ll2007 .l2008aUTriDathi et al. k200al and nume ricallv (e.g.. 



Gilbert et 



Manchester et al 



al. 



2006; Birn et al. 20061 : iFan Gibson l2007l 'l . especially in the con- 



text of eruptions driven by the kink instability. In particular, iFanI (j2005l ) demonstrated that the 
evolution of the kink instability facilitates the loss of confinement of a flux rope. Specifically, the 
writhing rotation changes the direction of the rope apex r elative to the overlyin g arcade, which 
helps the flux rope to "rupture" through the arcade field. I Gibson &: FanI (j2006l ) extended Fan's 
simulation in time to show that the flux rope bifurcates into two, with one part erupting and the 
other staying behind. In the latter simulation, one can see that overlying loops are pushed upward 
and aside as the flux rope kinks and expands, and that following the rupture of the arcade, those 
overlying loops that have not reconnected with the rope field quickly contract back to the core field, 
which is reasonable, considering the reduction of the magnetic pressure in the c ore field due to the 
rope e scaping. This is reminiscent of the "implosion" conjecture proposed by iHudsonI (j2000l ). to 
which iLiu et al.l (|2009al ) lent support in observation. In that case, however, the reduction of the 
magnetic pressure is presumably due to the release of the free magnetic energy during the early 
phase of the flare. 

In this Letter, we present observations of the coronal loop contraction in the wake of a filament 
eruption. The eruption occurred in the active region NOAA 9502 (S27E43) on 2001 June 15, 
associated with a GOE S class M6 flar e start ing at about 10:00 UT. By studying the photospheric 
helicity transportation, iRomano et al.l (j2005l ) concluded that the reverse S-shaped filament erupted 
due to the accumulation of magnetic helicity exceeding the kink instability threshold. Here we 
concentrate on a late-phase hard X-ray (HXR) burst recorded by Yohkoh (see the Yohkoh flare 
catalogu^, which occurred about 26 minutes later in the same active region, following the loop 
contraction. We will analyze the observations in Section 2, discuss the results in Section 3, and 
make concluding remarks in Section 4. 



2. Observation and Analysis 

2.1. Coronal Implosion 

The reverse S-shaped filament is obviously visible in Figure[T]^a) as a dark, absorption structure 
in the bright active region, as observed with the TRACE 195 A filter. A bundle of brightening 



" jhttp: //gedas22 . stelab.nagoya-u. ac . jp/HXT/catalogue/index.html[ HXT event number: 26710 
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Fig. 1. — Snapshots of the expansion and contraction of coronal loops associated with a filament 
eruption. Images are observed with the TRACE 195 A filter. The loops of interest are marked by 
white arrows. 
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a) TRACE 195 A Slit Images 
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b) GOES X-Ray Flux 
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d) Yohkoh Lightcurve 
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Fig. 2. — Coronal implosion observed in EUV and the lightcurves of the associated flare, (a) Slices 
of TRACE 195 A images cut by the slit as shown in the inset are placed on the time axis. The 
TRACE image in the inset has been rotated clockwise for 145°. Labels indicate speeds derived 
from the apparent "slopes", in km s^"*^. (b) COES soft X-ray fluxes in 1-8 A (grey) and 0.5-4 A 
(dark), (c) Time derivatives of COES fluxes, (d) HXT counts per second and per subcollimator 
(SC) are scaled by the y-axis on the left, and SXR counts recorded by the Soft X-ray Telescope 
(SXT) Be-119 filter (integrating over the post-flare arcade of the M6 flare) are scaled by the y-axis 
on the right. Dotted lines mark the time interval for the HXT image synthesis. 
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loops astride the filament are marked by an arrow. The loops are supposed to connect the opposite 
polarity separated by the filament channel (c.f. §2.2). From about 09:57 UT to 10:06 UT, as th e 



filament ascended slowly with a counterclockwise rotation of its axis (also see I Green et al.l 120071 ) , 
the overlying loops are observed to rise with it (Figure [II|a-c) ) , presumably being pushed upward 
by the filament, but to contract (Figure [T]|^d-f)) from 10:11 UT onward, after the filament "broke" 
through the arcade at about 10:08 UT (Figure [TJcy); a comprehensive animation of TRACE images 
for this event is available at the TRACE websitqll). Lightcurves of the flare associated with the 
filament eruption are presented in Figured! Within one hour of the eruption, one can see that the 
filament channel in Ha was filled with dark material again (Figure EU^d)), and that post-flare loops 
are shown as dark features perpendicular to the channel, suggestive of the partial eruption nature 
of the event. 

To study the loop evolution in detail, we rotate TRACE images clockwise by 145° so that the 
loops of interest are positioned upward (see the inset of Figure [2^a) ) . A slit of 290 pixels long and 
10 pixels wide is placed across the group of coronal loops, and the slices of TRACE images cut by 
the slit are placed on the time axis in Figure [2^a). Initially, the loops expanded slowly at about 
8 km s~^, and then at a faster speed of about 56 km s~^, roughly the same speed of the rising 
filament. At about 10:06 UT, the filament started to rise upward explosively, simultaneously with 
the onset of the flare impulsive phase, as indicated by the time derivatives of C OES soft X-ra y 



(SXR) fluxes (Figure [2h), which usually resemble the HXR profile (Neupert effect; lNeupertlll968l ). 
At 10:07:55 UT (Figure [^c)), the filament already lost its coherent shape, and the slit was filled 
by dark filament material during 10:06-10:08 UT (Figure [2]^a)). The group of loops became very 
diffuse and faint due to shorter exposures with enhanced brightening in the flaring region. But one 
can still see that as early as about 10:11 UT, the loops began to contract inward at various speeds 
ranging from 60 km s~^ to 140 km s^^. As an instruc tive comparison , the contraction of coronal 



loops observed in the flare early phase is at ~5 km s~^ ( Liu et al.ll2009al ). while the shrinkage speed 
of relaxing cusp-shaped field lines observed in SXRs is only -^2 km s~^ ( Linll2004l . and references 
therein). 

The contraction is observed to sustain for at least 12 minutes from 10:11 to 10:23 UT. During 
that period the loops of interest contracted from the maximum projected height of ^ 1.3 x 10^'^ 
cm, with respect to the flare footpoints (see Figure El^e)), to a height of ~ 6 x 10^ cm, where 
they became hardly visible, being overwhelmed by the bright active region. The post-flare arcade 
originating from the M6 flare is about 4 x 10^ cm high (projected). At an average speed of 100 km 
s~^, the contracting loops are expected to be stopped above the arcade in > 200 s, i.e., at about 
10:26 UT, coincident with the HXR burst. 
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Fig. 3. — Multi- wavelength observations of the M6 flare and the late- phase HXR burst, (a) A 
pre-flare Ha image taken by the the Kanzelhohe Solar Observatory (KSO) is overlaid with an line- 
of-sight magnetogram obtained by the Michelson Doppler Imager (MDI) on board SOHO at 09:35 
UT. Contours levels are 100, 200, 400 and 800 G for positive polarities (reddish orange), and -800, 
-400, -200, and -100 G for negative polarities (sky blue) (b) A Ha image shows the classical two- 
ribbon flare associated with the filament eruption, (c) HXT and SXT contours for the late-phase 
burst are overlaid on the KSO Ha image taken at approximately the same time. The contour levels 
are 10, 20, 40 and 80% of the maximum brightness, (d) shows a post-flare Ha image, (e) The same 
contours as in (c) are overlaid on a TRACE 195 A image taken at approximately the same time. 
The field of view is the same for panels (a)-(e). (f) shows a full-disk SOHO EIT 195 A image taken 
at 10:36 UT. A brightening at the east limb is marked by an arrow. 
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Table 1. Spectral Characteristics of the late-phase HXR burst 



Algorithm 


L/Ml M1/M2 


^SFP 

L/Ml M1/M2 


^NFP 

L/Ml M1/M2 


T^S (MK) 
L/Ml M1/M2 


MEM 
Pixon 


3.7±1.1 
3.1±0.7 


1.6±0.4 3.2±0.3 
1.8±0.3 2.6±0.2 


3.1±0.6 
2.3±0.3 


163±64 
218±94 



Note. — Spectral information from HXT images is derived by calculating flux ratios between 
adjacent energy bands. A region growing method is used to specify the integration region for 
each individual HXR source fsee iLiu et aljBoOQbl . for details). CS stands for coronal source. 
SFP (NFP) is referred to the footpoint source located in the south (north). 



ft 




Fig. 4. — Cartoon of a shrinking trap as a simplification of the observations. Top panel shows that 
initially the loop is pushed upward. Bottom panel shows that the loop shrinks toward the post-flare 
arcade. Lq is the loop's initial height (~ 8 x 10^ cm) measured from TRACE images. Similarly, 
Lmax is the loop's maximum height (~ 1.3 x 10^'^ cm), and I the height of the post-flare arcade 
(~ 4 X 10^ cm). V is the average contraction speed. Bi and denote the magnetic field at the 
height of / and Lq, respectively. The dotted line indicates the loop position at 10:23 UT. 



- 8 - 



2.2. Above-The-Loop-Top Hard X-Ray Source 

The HXR burst observed by Yohkoh HXT features a single peak that lasts for about 2 min. 
The four HXT bands show a similar profile (Figure EJd)), with no corresponding rise in SXRs 
(Figure [2jb) and (c)). The dotted lines in Figure [2]||d) mark the time interval of data accumulation 
for reconstructing HXR sources. The images synthesized with the Maximum Entropy Method 
(MEM) are shown as contours in Figure [3] (c) and (e) , which are very similar to those reconstructed 
with the Pixon method (therefore are not shown here). Note three minutes later (~10:29 UT), 
an M2 flare was recorded by GOES, which, however, is pinpointed to a different active region 
(NOAA 9506), in view of the following observations: 1) The SOHO FIT 195 A image taken at 
10:36 UT (Figure [3l^f)) shows a brightening point at the east limb (marked by an arrow), with no 
corresponding brightening in the previous FIT image at 10:24 UT (not shown); 2) HXT L- and 
Ml-band images (not shown) synthesized during the M2 flare also exhibit a point source co-spatial 
with the 195 A brightening; 3) The SXR counts integrating over the post-flare arcade of the M6 
flare decrease exponentially during the M2 flare (Figure [2^d)). 

The most striking feature of HXR burst is the existence of a coronal HXR source whose 
centroid position is about 15" above the SXT loop top (Figure EJe)). No significant SXR emission 
was detected at the coronal source location, indicating a rather low thermal plasma density. In 
addition, the coronal source is so lacking in L-band emission that no plausible spectral index can be 
derived from the L and Ml channel ratio (Table [T|). In fact, the whole HXR burst is also relatively 
lacking in L-band photons (see Figure [2]J^d)). This indicates a short of low-energy electrons, which 
normally would have constitute the bulk of the total energy. Thus, the energy input into the 
chromosphere may be too small to drive chromospheric evapo ration, which explains the absence of 



the Neupert effect (|Veronig et al.l 120051 : IWarmuth et al.l 120091 ). 



On the other hand, the spectral index of the coronal emission derived from the Ml and M2 
channel ratio is rather hard, similar to those of the footpoints (Table [T|). As a result, thermal in- 
terpretation cannot yield plausible temperatures. One may wonder whether the presumed coronal 
source is actually a third footpoint. However, this is not supported by the bipolar magnetic con- 
figuration involved in this classical two-ribbon flare (Figure [3|^a)-(c)), nor by the apparent bipolar 
coronal loops in the active region (Figured]). Moreover, one can see in Figure [3jc) that the conju- 
gate footpoints are located at the two flare ribbons, while no visible Ha brightening is observed at 
the coronal source location. 



If we assume the thin-target hypothesis for the coronal source {6'^^ = 72"^ — l|||), and thick- 
target for the footpoint emission {6^^ = + 1), then the injecting electron population for 



^http: //trace . Imsal . com/POD/bigmovies/f ilainents/T195_2001Junl5New.mov | 

■''We denote the photon spectral index derived from the L- and Ml-band ratio as 71, and that from the Ml- and 
M2-band as 72, with the superscripts, CS and FP, standing for the coronal source and the footpoint, respectively. 
Corresponding electron power-law indices are denoted by 5 in a usual fashion. 
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the coronal source, S^^ , is hardening by about 1.5 powers, relative to that for the footpoint, 
5^^, which is consistent w ith electron trapping in the weak diffusion limit (c.f., the Apendix in 



Metcalf k Alexandeilll999l l. 



3. Discussion 



The characteristics of t he above-the-loop-t op source presented in this Letter are reminiscent 
of the fa mous Masuda flare (iMasuda et al.lll994l ). which is regarded as a very unique eve nt (see the 
review bv lKrucker et al.ll2008l ). In a recent revisi t of the Masuda flare ( Liu et alj|2009bl ). however, 
we noticed that in the calibrated HXT images (jSato et al.lll999l ). the coronal source looks like a 
normal loop-top source, located slightly above the SXT loop top but within the 10% brightness 
contour of the HXT L-band th ermal loop, rather than d istinctly above the loop top as in the original 
HXT images (e.g., Figure 3 in lAschwanden et al.l 119961). Accordingly, its spectral indices also fall 



Tomczak 



2001 



Krucker &: Linll2008l ) . These changes 



into the normal range of coronal emission (4-7; 
have relieved, in large measure, the longstanding difficulties of understanding the original Masuda 
flare, viz., the lack of L-band emission of the coronal source, and the low density at the coronal 
source location. 

On the other hand, the so called above-the-loop-top sources reported in the R HESSI era have 
been primarily, if no t exclusively, found in a double coronal source morphologv (ISui Holman 



2003 



Sui et al 



2004 



Pick et al.l l2005l : IVeronig et al. 



2006 



Li &: Gan 



2007 



Liu et al 



2008b!). In 



that case, the lower coronal source is located at the thermal loop top, while the upper source can 
be located as far as 3 0^^ above the lower one, moving upward at a speed as high as 300 km s~^ (e.g., 
Sui &: HolmanI I2OO3I ) . The two sources "mirror" each other with respect to a presumed X-point 
reconnection site, in that higher energy emission com es from lower altitudes for the upper s ource, 



while the lower source exhibits a reversed order (e.g., ISui &: Holman 



spectra of the two sources have a similar power-law index (|Liu et al 



2003 



2008bl ). suggesting that the 



Lhi et al.ll2008bl ). The 



injecting populations are intimately related, presumably accelerated at the X-point. If an X-point 
configuration were assumed below the above-the-loop-top source studied here, as illustrated in 
Figure [3]^c), an additional loop-top source should also be expected, which is obviously not the case. 

Alternatively, the contracting loops may be relevant in the HXR production, since a shortening 
of the field line length indicates a reduction of the magnetic energy intensity, which is available as 
free energy to heat and accelerate particles. This effect has been discussed in the context of newly 
recon nected cusp-shaped lines relaxing into potential-like lines at lower altitudes (jSomov &: Kosugi 
1993). 



As a simplification of our observation. Figure [4] illustrates that a magnetic loop is initially 
stretched and pushed upward (top panel), and later forms a shrinking trap (bottom panel) in 
which mirroring particles are accelerated in a way similar to the trap-without-shock scenario in 
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Somov &: Kosugil (|l997l ). The lifetime of our shrinking trap, 

^ Lmax - I 9 X 10^ cm 



107 



900 



cm s 



where v, the shrinking speed of the coronal loop, is approximated to be 100 km s ^ (see Figure [2]^ a)). 
The collisional deflection time of an electron of energy E is given by. 



t{E) = 0.95 



1 keV J 



20 
huV 



,10^ cm~3 

where the Coulomb logarithm, InA = ln(8 x 10^ T^Ue ^^'^) — 20 under typical coronal conditions. 
Thus, for 10-100 keV electrons in the coronal trap with a density in the order of 10^ cm~^, the 
trapping time is of order 100-1000 s, comparable to the trap lifetime. 

For an 1-10 keV electron, its thermal speed — (2 — 6) x 10^ cm s~^, so its characteristic 
time between two subsequent reflections, r < 2Lmax/ve — (4 — 13) s. Since t <^T, the conditions 
of the periodic longi tudinal moti ons change adiabatically. The longitudinal adiabatic invariant is 
therefore conserved (|Somovl 120061 . Chapter 6), i.e.. 



P\\{t)L{t) = const. 

Since L^ax/l — 3.3, when t ^ T, P\\{t) 3.3P||(0). In addition, the transversal adiabatic invariant 
is conserved in the non-relativistic approximation. 



Plit) 
B 



const. 



Assuming that both the initial and the final state of the magnetic loop are qu asi-potential, we 
estima te the field strength at the loop top, following the empirical formula from lDulk &: McLean 
hdrd ). i.e., 

Bi ^0.5 {I /Rq)-^^^ ^36 G, 

and Bl ~ 0.5 (Lq/Rq)'^^"^ ~ 13 G, where Bi and B^ denote the magnetic field at the height of 
/ and Lq, respectively. As the contracting loop approaches the post-fiare arcade (bottom panel 
of Figure HD, B Bi, hence P±{t) 1.7P±{0), and Ekit) = {P±itf + P||(t)2)/2me ^ 14Sfc(0). 
Thus, hot electrons of several keV can be accelerated to nonthermal energies in this trap. 



4. Conclusion 

We suggest that the escaping of the kinking filament in the 2001 June 15 event results in the 



( Hudson 


2000: 


Liu et al. 


2009a ^ 



efficiently accelerate electrons. This explains the HXR burst during the flare decay phase, which 



-Il- 



ls characterized by nonthermal coronal emission in consistency with electron trapping in the weak 
diffusion limit. Although magnetic reconnection cannot be excluded in this picture, the shrinking 
trap model seems to offer a more self-consistent interpretation. 

We expect this acceleration process to occur generally in eruptions in which the arcade field is 
only partially opened. Free energy builds up in the arcade when field lines are pushed upward by 
the eruptive flux rope, since the magnetic energy associated with the loop. 



Svr 



B^ds dl, 



where s is the loop cross section, assuming to be constant, and / the loop length. As the stretched 
field relaxes in the wake of the rope escaping, trapped particles gain energy from the increase in par- 
allel momentum, due to shortening field lines, and from the increase in perpendicular momentum, 
due to strengthening magnetic field at lower altitudes. 

This process may play a role in some of the HXR bursts observed during the flare late 
phase, which often show weak soft emission, hardening spectra, and/or high co ronal sources (e.g., 
Frost fc Denni] Il97ll : iHudsonI llQTsI : lOrigis fc Ben3 boosl : Iwarmuth et^bood ^. Such late-phase 
bursts have been o ften attributed to the acceleration and trapping of electrons in the post-flare 



loop systems (e.g.. ICliver et al.lll986l ). or, to the acceleration of electrons in a shock front (e.g., 



Frost &: Dennis! Il97ll ). The mechanism suggested here, however, involves a transfer of the free 
magnetic energy from the core field to the confining arcade during the eruptive process. The en- 
ergy transferred is then made available to the trapped particles in the aftermath of the partial 
eruption, through the contraction of the loops that have not reconnected. 
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